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ABST RAGA 


The GFCS MK86 solves the ballistics ShoRiee for the 5'/54 
gun iteratively utilizing table look-ups. Each ballistics parameter 
table is sub-divided two-dimensionally, corresponding to selected 
elevation and slant range intervals. The final value for each para- 
meter is found by interpolating between parameter values evaluated 
at elevation angles which bracket the actual elevation angle. 

This thesis describes the method of determining the ballistics 
solution for AA cartesian ballistics. A description of the ballistics 
geometry and of the model utilized is presented, together with a 
FORTRAN simulation of the MK86 computer program and a listing 


of the ballistics tables used in this system. 
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I. INTRODUCTION 


A. GENERAL DESCRIPTION 

The GFCS MK86 is the first digital gun-fire control system 
to be considered by the U. S. Navy. It is designed to control the 
5'/54 gun-mount and to solve the fire control problem for the air, 
surface, shore-direct or shore-indirect modes of operation. In this 


thesis only the AA cartesian ballistics mode is considered. 


5, PROBLEM DESCRIPTION 

A computer simulation program was felt to be a prerequisite 
to a complete study and analysis of the GFCS MK86 AA model This 
in turn necessitated a careful study of the existing documentation 
and machine-language software of the MK86 system. 

The FORTRAN simulation of the MK 86 AA mode was divided 
into three sections: The pre-filter processing section, programmed 
mest J. D. GORMAN; The three-dimensional numerical filter, by 
iT T. J WILL; and the ballistics section contained in this thesis. 
These three programs can process radar target information from 
target detection to ballistics output prior to the generation of gun 
orders. 

Another objective of this project was the generation of adequate 
documentation of the system. This phase was also divided into the 


three sections described in the previous paragraph. 
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(oe SUMMARY OF INTRODUCTION 

This thesis contains a description of the MK86 ballistics program 
including the ballistics intercept geometry, ballistics model used, and 
sequence of ballistics computations. A description of each subprogram 
indicating its function is also presented. The FORTRAN simulation 
program including the ballistics tables is listed and differences 
between the machine-language program and the FORTRAN program 


are indicated. 


Il. INTRODUCTION TO MK86 BALLISTICS 


The GFCS MK86 ballistics program accepts stabilized and 
filtered target position and rate data, together with ownship's data, 
and computes the gun ballistics required for the generation of gun 


pointing and fuze orders. 


As. BALLIS FICS SETUP (1.%. INITIALIZATION) 

The ballistics parameters suchas air temperature, air density, 
wind, and initial velocity are entered at the console and corresponding 
increments from noninal values are computed. 

Target and ownship's data are used to compute the target's 
present position and rate relative to the gun's position at ballistics 
time. This position of the target is called point 1 (see Figure l), 
and is computed in spherical coordinates. The relative rates are 
resolved into range and crossrange components. 

A test is made to determine whether the input data is associated 
with a new target. If it is, an acquisition mode is declared and the 
processing sequence calls the acquisition subprogram. Inthe case 
of an old target, track mode is designated and the acquisition sub- 
program is bypassed. This is illustrated in the flow chart of 


Figure 2. 
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Bt BALLISTICS ACQUISITION 

Before a ballistics solution can be developed, a first estimate 
of future target position is required. This is called the advance 
position, or point 3 (see Figure 1). If the system has been tracking 
the target, then the last ballistics solution, updated, serves as the 
advance position. With a new target, however, some other means 
for estimating a reasonably accurate point 3 is necessary. This is 
the function of the ballistics acquisition subprogram (see Figure 2). 

In the acquisition phase, point 3 is computed as follows. The 
projectile time of flight to the present target position (point 1!) is 
computed from the ballistics table in the normal way (see Section III). 
Using this time of flight and the target's present position and velocity, 
an estimate is made of the target's future position (point 3). A number 
of logical tests are employed to check the validity of the estimate. 
The slant range, bearing, and elevation of point 3 are designated 


io, bY3, and 3 respectively. 


C. DETERMINATION OF POINT 2 

The ballistics program next uses the advance position of the 
target (point 3) as the ares point and computes a predicted "impact"! 
point from the ballistics tables. This is called point 2 and the 
corresponding slant range, bearing, elevation, and time of flight are 


designated RZ, BY2, E2, and T2, respectively. These values differ 
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from the point 3 quantities due to non-nominal powder and air tem- 


peratures, and because of coriolis, drift, and wind effects. 


ioe CALCULATION OF POINTS 21 AND 22. 

The target position which corresponds to the computed time of 
flight T2 is calculated and designated as point 21 (see Figure 1). Then 
the differences in slant range, bearing and elevation from those at 
point-c (the predicted projectile impact point) are computed. These 
increments are designated IR21, IBY21, and IE2Z1, respectively, and 
a corresponding linearized estimate of the time of flight increment 
IT2 is computed. 

Finally the point 21 target position is updated by the time 
increment ITZ to determine point 22 (see Figure 1), which is the 
first estimate of the point where the target and projectile will coincide 


in both position and time. 


i. BALLISTICS ITERATION (I.E. SECOND PASS) 

After point 22 noe been calculated once, as indicated in the 
ee sections, the program uses point 22 (modified) as a new 
(improved) point 3 and goes through the calculation to point 22a 
second time. This is shown in Figure 2. Finally the increments 
imncz, IBY22, 1E22, and ITZ2 between point Z and point 2Z are 


computed for use in the ballistics output program. 


12 


i. CLOSURE TESTS 

A test is made at the end of each of the two passes described 
above to determine if the final point 22 is a valid (i.e. accurate) 
solution. This is done by comparing point 21 data and point 2 data, 
which is called the closure test. The test after the first pass merely 
sets a flag, which is used to modify the criterion used in the test 
which follows the second pass. If the second-pass solution fails to 


close, then the computation aborts for this time interval. 


1s BALLISTICS OUTPUT 

If closure has occurred, the gun elevation, bearing and time of 
flight which were calculated with point 2 as the aim point are modified 
by the amount necessary to move the impact point to point 22. This 


data is then transmitted to the gun and fuze order programs. 
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Figure 2 
FUNCTIONAL FLOW CHART 
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WI, BALLISTICS EFFECTS MODELING 


The MK86 ballistics modeling program provides the means for 
determining the values and effects of the various ballistic parameters. 
The approach taken and the structure of the modeling is described 
here for AA cartesian ballistics (see Figure 3). The basic function 
of the model is to provide projectile impact or intercept point parameters 
based on advance target position parameters (that is, to provide point 2 
output data from point 3 input data). The ballistics parameters are 


broken down into the three categories described below. 


2s CATEGORY-ONE PARAMETERS 
be Gun elevation angle, EG 


Die Hime of flight, V3 


BY Range effect per unit initial velocity variation, DRDU 

4. Range effect per unit air density variation, DRDH 

5. Range effect per unit air temperature variation, DRDTF 

6. Time of flight effect per unit initial velocity variation, 
DT3DU 


ie Time of flight. effect per unit-air density variation, DT3DH 
Gc. Time of flight effect per unit air temperature variation, 


D3 DLE 


Ig 


thd Stiel 


Each of these parameters is modeled by a cubic polynominal 
function of range of the form: 
F=K0+K1(IR3)+K2(IR3)° +K3(IR3)> 
=K0+IR3(K1+IR3(K2+K3(1IR3))) 
Where F is the parameter to be determined and the K's are coefficients 
obtained from tables which are arranged in elevation ''stacks", each 
of which is 'partitioned'’ according to slant range. This is illustrated 
by the table below: 


TABLE, I 
Typical Ballistics Table for One 


Parameter (such as Time of Flight) 


Stack 1 Stack 2 Stack N 
Stack Stack Same format 
elevation elevation as other 
angle #1 angle #2 stacks 
List of List of 
range range 
intervals intervals 
in this in this 
stack stack 
Coefficients Coefficients 
for first | fom fis t 
range interval range interval 
Coefficients Coefficients 
for 2nd range for 2nd range 
interval interval 
etc. etc. 


The quantity IR 3 in the above cubic polynominal is the incremental 
slant range (i.e. the slant range of interest minus the lower range 


bound of the range partition in use). 
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That is, each category-one parameters is modeled by a cubic 
equation in each range-elevation partition, as shown in Figure 4. It 
iepthe coefficients of the cubices oe are provided by the ballistic 
tables. Inthe simplest case, two values of F are computed, one 
corresponding to each of the elevation stacks which bracket the actual 


elevation. The final value of F then follows from linear interpolation. 


B. CATEGORY. TWO PARAMETERS 

is Coriolis effects on range, crossrange and time of flight 

(ie Drift effect 

ae Angle of fall factor and calculation of tangent of the angle 

of fall 

Each of these parameters is modeled by an analytic expression 
modified by a correction factor which is calculated in essentially the 
same way as a category-one parameter. In this case the factor K is 


computed by a linear incremental slant range relation, ‘K=C0+C1(IR3). 


. CATEGORY-THREE PARAMETERS 


Ie Range effect per unit projectile weight variation 
2. Time of flight effect per unit projectile weight variation 
oN Wind and motion effects on range and crossrange 


These parameters are calculated in special subroutines which 


do not require range and/or elevation dependent table lookup. 
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Figure 3 
BALLISTIC MODEL 
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Figure 4 
PARAMETER AS A FUNCTION OF ELEVATION AND RANGE 
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IVa (RHE BALLISTICS: TABLES 


The ballistics tables were derived by modeling range table 
data provided LEC by NWL Dahlgren [Ref. 1]. There are six sets of 
tables. Each table corresponds to one or more ballistics parameters. 
The AA tables are of the following form: 
by. Each parameter's section of the table is divided into 
elevation stacks. 
ae Each elevation stack is divided into range partitions 
or intervals. 
De For each range partition there are a set of coefficients 


which are used to evaluate the chosen parameter. 


A, MK86 MACHINE LANGUAGE TABLES 

The tables are located in the computer in sequential order. The 
address of the first entry for each parameter is stored in the ETABAD 
matrix. This address is.also stored in ISRCH(1). To find the address 
of the first elevation stack to be used, the address stored in ISRCH(1) 
is modified by the contents of the PUSHTB matrix. The first 8 bits 
of each PUSHTB entry contains the number of entries to be bypassed 
in order to locate the address of the initial elevation stack for the 
low elevation AA mode (advance elevation less than 35°). The last 
8 bits indicate the number of entries to be bypassed for the high ele- 


vation AA mode. The modified address is then stored in ISRCH(Z2). 


ZS 


The first entry of each elevation stack contains: 


be Stack elevation angle 
Ze Number of range partitions, minus one 
Si. Internal shift count 


The stack elevation angle is stored in STACK(3), the number of range 
partitions is stored in the NP matrix, and the internal shift count 
necessary to position the binary point in the KSH matrix. 


A complete listing of the tables may be found in reference 2. 


B. BALLISTICS TABLES IN THE FORTRAN SIMULATION 

The entry into the the organization of the tables is different in 
the FORTRAN simulation program from that in the machine language 
program. 

The ballistics tables are stored in a three-dimensional array. 
The coordinates of he matrix are parameter, elevation stack angle, 
and the location of a particular item inthe stack. The first entry in 
a particular elevation stack is the elevation stack angle, and it is 
followed by the number of range partitions. Then the limits of the 
range partitions are listed, followed by the coefficients for each 
range partition. The location ina stack of the first coefficient 
required is calculated by multiplying the number of coefficients times 
the number of the range partition of interest and adding the location 


of the first range partition plus 3. Then the location is stored in T17. 


A listing of the ballistics tables starts on page 90. In this listing, 
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each table corresponds to a ballistic parameter (see Table 2), and 


each column corresponds to an elevation stack. 


Table 


10 
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TABLE 2 


Ballistic Tables (FORTRAN) 


Parameter 

Time of flight (T3) 
Super-elevation angle (EGB) 
Range effect due to change in 
initial velocity (DRDU) 

Time of flight effects due’to 
change in initial velocity (DT3DU) 
Range effect due to change in air 
density (DRDH) 

Time of flight effects due to 
change in air density (DT3DH) 
Range effect due to change in 

air temperature (DRDTF) 

Time of flight effects due to change 
in air temperature (DT3DTF) 
Tangent of angle of fall (TANW) 
Coriolis effects 


Drift effects 
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V. BALLISTICS SUBPROGRAMS 


Since this project was concerned only with the AA aspect of the 
GFCS MK86 system, all references to surface ballistics, grid mode, 
ROS mode, TDT mode, and spherical ballistics have been omitted. 
For the flow chart of the subprograms see Figure 5. A complete 


explanation and flow chart may be found in References 3 and 4. 


A. TARGI 

This subprogram freezes or holds the console inputs for the 
duration of the % second ballistics computation cycle, and calculates 
and sets the following indicators: 

i, AA ballistics 

as Cartesian ballistics 


She Ammo type, charge 


4. Gun and gun control console assignments 
5. Target assignment 
6. Acquisition or tracking ballistics 


In addition, TARG1 computes non-nominal increments that may 
exist in various ballistics parameters, performs a ship's course 
update, and determines gun parallax corrections. Calculations are 
also performed to update the target data to "ballistic time" (i.e., the 


next 'integer''% second). 
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B. TARG2 

This subprogram performs the acquisition ballistics computations 
for AA cartesian ballistics. Its function is to calculate point 3 (the 
theoretical aiming point) from point 1 (the target position at ballistics 
time). This subprogram determines the first advance position (aim 
point) when no valid prior advance position exists. This situation 
would occur whenever: 

i A new AA target has been assigned 

2. There was no ballistics solution during the previous 


*%-second ballistics computation cycle 


Si There has been a change in the mode assignment for 
the gun 
4. There has been a change in the ammo type.or charge 


The countdown flag (KGOOD) is set to one, indicating that the 
acquisition subprogram must be completed more than once before the 
gun is released from any previous state it may have beenin. The 
meguire flag is set a one to indicate the acquisition ballistics mode. 

The advanced elevation (E3) is set equal to the elevation angle 
at point 1 (El). TMONIT, which controls the ballistic table computa- 
tions, is then entered to find the nominal time of flight (Tl) from RC, 


tie*range to point I. 
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8 BCOMPA 

This subprogram computes the time of flight to point 3 (T3) 
and its derivatives with respect to range and elevation. It also 
computes the mixed second partial derivative with respect to range 
and elevation. Using E3 and R3 as inputs, TMONIT is called and 
the time of flight is returned in TEMP2Z, the range derivative (DT3DR) 
in TMP3, the elevation derivative (DE3DE) in DPDEP, and the mixed 


partial (DT3RDE)in DPDED. 


Dy, BCOMPB 

This subprogram computes the super elevation angle (EGB) and 
its derivatives with respect to range and elevation, as wellas the 
mixed second partial with respect to range and elevation. TMONIT 
is called and the super elevation angle is returned in TEMP2, the 
range derivative in TMP3, the elevation derivative (DEGDE) in 
DPDEP, and the mixed partial (DEGRDE) in DPDED. Then EGB and 
E3 are summed to form the total elevation angle E4=E3+EGB, which 


is limited to a maximum of 85 degrees. 


BE. BCOMPC 

This subprogram computes the derivative of range with respect 
to a non-nominal increment of projectile velocity (DRDU) and the 
derivative of DRDU with respect to elevation (DRDUDE). Then the 
actual range effect (URH) due to non-nominal projectile velocity is 


computed. 
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iy. BCOMPD 

This subprogram computes the derivative of the time of flight 
with respect to the increment in projectile velocity (DR3DU) and the 
time of flight effect (UT3) due to non-nominal projectile velocity. 
The mixed second partials of time of flight with respect to increment 
in projectile velocity and elevation (DT3UDE) and range (DT3UDR) 
are also computed. 
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G. BCOMPE 

This subprogram computes the derivative of range with respect 
to non-nominal increments in air density (DRDH) and the derivative 
of DRDH with respect to elevation (DRDHDE). Then the actual range 


effect due to non-nominal ballistic air density (HRH) is computed. 


i. BCOMPF 

This subprogram computes the derivative of time of flight with 
respect to a non-nominal increment in air density (DT3DH) and the 
actual time of flight effect (HT3) due to DT3DH. The mixed second 
Peeiats of time of flight with respect to air density and range and 


elevation (DT3HDR, DT3HDE) are also computed. 


i? BCOMPG 
This subprogram computes the derivative of the range with 


respect to a non-nominal increment in surface air temperature (DRDTF), 
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the actual range effect (TFR3), and the mixed partial of range with 


respect to temperature and elevation (DRTFDE). 


e BCOMPH 

This subprogram computes the derivative of time of flight with 
Bespect to a non-nominal increment in surface temperature (DT3DTF); 
computes the actual time of flight effect (IT FT3) due to air temperature; 
and calls the weight effect subprogram (RTWT). 

The subprogram then performs a partial sum of all the range 
effects and time of flight effects calculated so far and calculates the 
derivative of the time of flight to point 2 (the predicted impact point) 


with respect to range (DT2ZDR). 


K. BCOMPI 

This subprogram computes the downrange and crossrange 
components of gun motion (DMR30G DMGOQG) before calling GUNMOT, 
the gun motion subprogram. It then computes the tangent of the angle 


of fall (TAN W) and calls the wind effects subprogram, WIND. 


teed BCOMPI 
This subprogram obtains a first partial sum of range effects, a 
second partial sum of time of flight effects, and calls CORLIS to 


obtain coriolis effects. 
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M. BCOMPK 

This subprogram computes the intercept geometry by determining 
the horizontal range, vertical range, slant range, and bearing to the 
predicted projectile impact point (point 2). It also computes the 


elevation angles for the gun. 


N. TARG3 

This subprogram performs the iteration and closure computations. 
Its function is to update through two iterations the previous advanced 
aim point (point 3) such that satisfactory coincidence is obtained 
between the predicted impact point (point 2) and the predicted target 
position at impact time (point 21). The closure equations determine 
when there is sufficient closeness between the two points for an 
acceptable solution. 

The input o8 Abe. to TARG3 come from two sources. In the 
tracking mode, the data from the previous pass through TARG3 to 
BCOMPK is used to update current values. In the acquisition mode, 
the values calculated in antes eae TARG2 through BCOMPK are 


used. 


oO. TARG4 
This subprogram computes the ballistics outputs. It computes 
the final inertial gun orders, gun order rates, and fuse time rate as 


well as the inertial impact coordinates and rates for eventual display 


ol 


(ty ° teryiTte» freia-ti of. In , 


on the B-scan. These computations are oriented to the start of the 


next 1/16 second time-slice. 


P, TMONIT 

For a given parameter, this subprogram does the initial setup 
and locates the proper table and elevation stack corresponding to the 
target's elevation. Then it calls appropriate subprograms to cal- 
culate the value of the parameter from the chosen table. 

In more detail, either two or three stacks are selected so that 
Reiz elevation values bracket the value of E3. Two are selected for 
parameters which are linearly interpolated with respect to elevation 
and three for those which are parabolicaly interpolated. The partition 
selection, coefficient look-up, and polynominal evaluation with respect 
to range are performed on each of the selected sets of coefficients 
yielding a value of the parameter and its derivatives corresponding 
to each of the bracketing elevations. The final value of the para- 
meter and its derivatives are then obtained by interpolating with 


respect to elevation. 


Q. AASRCH 

This subprogram controls the selection of coefficients for 
bracketing elevation stacks, the cubic polynominal evaluation, and 
the interpolation between elevation stacks. It also controls the 


calculation of the appropriate derivatives. For standard tables (four 


Sz 


coefficients per partition) the address ee the first coefficient is found 

by calling the SEARCH subprogram. Then the polynomial subprogram 
is called to evaluate the Ascites, which is stored, with its derivative, 
in the STACK matrix. Whenall parameters and derivatives correspond- 
ing to the selected elevation stacks have been evaluated, then the inter- 
polation subprogram is called to get the final values. For non-standard 
tables, only the address of the first coefficient is found for each selected 


elevation stack. These values are stored in the HOLDB3 matrix. 


Se SEARCH 
This subprogram, which is called by AASRCH, finds the correct 
range partition and the location of the first coefficient for that range 


partition. 


5. POLYN 
This subprogram, which is called by AASRCH, evaluates the 


polynomial from the coefficients and finds the derivatives with 


respect to range for each elevation stack chosen. 


oF INT ERP 

This subprogram performs linear or parabolic interpolation on 
a set of two or three parameter values with respect to elevation. 
The interpolates value is the value of the parameter at the actual 


advance elevation. 


DECSI 
This subprogram decodes the elevation information from each 


stack and stores it in STACK(3). 


‘eo SLIDE 
This subprogram moves the parameters stored in the STACK 


matrix up one location. 


WwW. FINDAL 
This subprogram finds the location of the next available elevation 


stack. 


X. CORLIS 

This subprogram computes ballistic Coriolis effects for slant 
range, altitude, time of flight, crossrange and downrange effects. 
The Coriolis table is a non-standard table with six coefficients per 
range partition. Three linear coefficients are computed and used to 


determine the effects of Coriolis. 


Ge DRIFT 

This subprogram computes the crossrange effect of projectile 
drift. The drift effect table is a non-standard table with two co- 
efficients per range partition. Two polynomials are evaluated for 
AA ballistics, one for each elevation stack bracketing the actual 


advance elevation. 
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Bis GUNMOT 

This subprogram calculates the ballistic effect due to gun 
motion (i.e. the advance range, time of flight, deflection and cross- 
range effects). Under conditions of no gun motion, the projectile 
leaves the muzzle with an initial velocity vector (U) at inertial true 
bearing BGY, and inertial elevation, EG. With gun motion, the 


component along the U vector will increase the magnitude of U. 


wee RE WT 
This subprogram calculates the (minor) effect of projectile 

weight on time of flight and range. Due to inertia, an increase of 
projectile weight will produce a decrease of initial velocity. It will 
also produce a decrease in drag acceleration, which is equivalent to 
a decrease in air density. Thus the range and time of flight effects 
due to a projectile weight variation can be computed as the sum of 
the equivalent effects of increments in initial velocity and ballistics 


air density. 


AB. WIND 
This subprogram calculates the total effects of wind on advance 


range, time of flight, and crossrange. 


AC. ZULIMIT 
This subprogram limits the Z coordinate of a target toa 
minimum value, based upon target horizontal range, earth's curvature, 


and height of the gun above the sea surface. 
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VI. FORTRAN SIMULATION 


Since only the AA mode was to be considered, all flags and 
variables necessary to indicate the AA mode were pre-set to insure 


that only the AA cartesian ballistics computations were performed. 


es ASSUMPTIONS MADE 
The FORTRAN program has no real time clock, therefore 
certain timing assumptions had to be made. These assumptions are: 
i, JOELCK (target valid data-time for ballistics) 


equals -1/32 second 


Bi, TBR (saved value of ballistic time) was initialized 
to 0 
Ds TB (ballistics time to next second) initialized to 


% second 
4. TBFIN (final ballistics time for next 1/16 second 


initialized to 15/32 second 


B. FORTRAN PROGRAM 
The program, which is listed starting on page closely follows 
the arrangement of the MK86 machine language program. It has 
been suitably annotated with comment cards, and otherwise should 
be self-explanatory. It has been written so that all ranges are in 


yards, and all angles are in degrees. 
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Three internal subprograms (DECSI, SLIDE, AND FINDAL) 
of TMONIT in the machine language program have been treated as 
separate subprograms in the FORTRAN program. Also, in BCOMPB 
the super-elevation angle is computed in minutes and then converted 


forge prees. 


Cc. VERIFICATION 

Complete verification of the FORTRAN program was not done 
due to lack of sufficient data. The separate subprograms were 
verified by calculations done by hand. Using data furnished by LEC, 
a partial verification was attempted. However, effects such as wind, 
coriolis, non-nominal parameters and ship's motion had been set to 


zero in the data provided [Ref. 5]. 
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